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Abstract 

A new rapid energization process within a supernova shock transition region (STR) is reported 
by utihzing numerical simulation. Although the scale of a STR as a main dissipation region is only 
several hundreds of thousands km, several interesting structures are found relating to generation of 
a root of the energetic particles. The nonlinear evolution of plasma instabilities lead to a dynamical 
change in the ion phase space distribution which associates with change of the field properties. As a 
result, different types of large-amplitude field structures appear. One is the leading wave packet and 
another is a series of magnetic solitary humps. Each field structure has a microscopic scale (~ the 
ion inertia length). Through the multiple nonlinear scattering between these large-amplitude field 
structures, electrons are accelerated directly. Within a STR, quick thermalization realizes energy 
equipartition between the ion and electron, hot electrons play an important role in keeping these 
large-amplitude field structures on the ion-acoustic mode. The hot electron shows non-Maxwellian 
distribution and could be the seed of further non-thermal population. The "shock system", where 
fresh incoming and refiected ions are supplied constantly, play an essential role in our result. With 
a perpendicular shock geometry, the maximum energy of the electron is estimated by equating a 
width of the STR to a length of the Larmor radius of the energetic electron. Under some realistic 
condition of Ma = 170 and Upe/^ce = 120, maximum energy is estimated to ~ 10 MeV at an 
instant only within the STR. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



Recently, the origin of the energetic electron related to astrophysical shocks has been 
investigated more and more from the point of view of the plasma kinetic processes [l|, 3, Is], 
y, y, y, [7, y, yi. within the theory of well-established classical diffusive shock acceleration 
(DSA) mechanism 10|, lUl] , a shock wave tends to be treated as a simple discontinuity. At a 
shock wave, however, we can find that many kind of plasma kinetics play a key role in the 
energy dissipation and generation of energetic particles. Correct treatment of the plasma 
kinetic process around the shock may succeed to settle so-called "injection problem" to DSA 
process (e.g. Ref.[l2|,Q,Q) as well as to capture a whole stage of DSA pro cess through 
nonlinear interaction between the shock structure and energetic particles 15|, ll6|, Il7|, Il8 |. 
Furthermore, not only as an assistant role to DSA mechanism, investigation of the energy 
release mechanism in terms of the plasma kinetics can be a powerful tool to search significant 
electron energization process in much smaller spatial/time scales than that DSA mechanism 
may need. 

In this paper, we focus on a nonrelativistic ion-electron supernova shock with the perpen- 
dicular geometry where the shock normal is perpendicular to the direction of the background 
magnetic field (e.g. Ref. 19|). Particle-in-cell (PIC) simulation has been carried out to inves- 



tigate coupling dynamics among the ion, electron, and field in the shock transition region. 
The shock transition region (STR) is the front side of the shock wave corresponding to the 
region fiow speed decreases gradually from the upstream toward the shocked downstream 
region. At the high Mach number shocks, part of the incoming ions are always refiected 
at the shock front and move into the upstream region. In detail, the STR means here 
the region from the leading edge of the refiected ion to around the end point of the first 
macroscopic ion gyration in a shocked region. The upstream electrons are immediately de- 
celerated just after they meet with the refiected ion. The velocity difference between the 
decelerated electron and the incoming ion is large enough to excite a series of micro-scale po- 
tentials through a strong two-stream instability which grows into a packet of large-amplitude 
spatial-oscillating magnetic field at the leading edge of the refiected ion (hereafter we call 
this "leading packet"). When the incoming ion begins to intermingle with the refiected ion, 
two ion components together make complex vortices in the velocity phase space, the electric 
and magnetic field structures change immediately. As a result, just after this dynamical 
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change, some components of the leading packet are converted into magnetic sohtary humps 
almost standing in the electron bulk flow frame. The electrons bouncing around between the 
leading packet and the standing magnetic humps get energy rapidly by the Ist-order Fermi 
type process. We define here, the Ist-order Fermi type acceleration as the acceleration of 
particles by a shock with an extended transition layer, while the original Ist-order Fermi 
mechanism (DSA) assumes an infinitesimally thin transition layer. Among magnetic humps, 
the electrons are also bouncing around and sometimes getting energy in a stochastic manner 
which depends on mainly local motional electric field caused by a combination of magnetic 
field of the humps and local bulk flow variation due to the macroscopic ion gyro-motion. 
Within the STR, strong thcrmalization occurs to attain energy cquipartition between the 
electron and ion. The hot electrons contribute to keep above large-amplitude structures 
on the ion-acoustic mode. The electron energization mechanisms stated above are quite 
dynamical process characteristic of a " shock system" . 

In the next section we introduce the simulation setup. Section III. presents simulation 
results: A. macroscopic overview of shock wave properties, B. the energy re-distribution 
process in the STR, and C. the structure of the STR and characteristic field generation. 
In subsection III.D we show electron energy spectra and discuss some examples of electron 
trajectories in the nonlinear field evolution. Section IV. summarizes and discusses our results. 

II. NUMERICAL SIMULATION SETUP 

One-dimensional (ID) electromagnetic relativistic particle-in-cell (PIC) simulation is 
adopted in the present paper. The shock wave is coUisionless and formed by so-called 
piston method. A high-speed plasma (velocity Uq) consisting of the electron and ion (its 
absolute charge value is equal to the electron) is injected from the left boundary (x = 0) 
and travels toward a; > region. Initially, a uniform magnetic field (strength Bq) per- 
pendicular to the flow direction (= z) is carried by the plasma. The y direction is deflned by 
the relation of each unit vector, x x y — z. At the right boundary, the plasma flow and all 
waves are reflected and accumulated to form a shock wave propagating toward the left. The 
right boundary of the simulation box is taken with enough distance from the left boundary 
to allow free shock propagation without any artificial downstream disturbances. The ini- 
tial plasma parameters are as follows: uq/c — 0.096, electron and ion beta (3e — Pi — 0.1 
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{f3j = SnTj/BQ, where Tj is the temperature for j species, j = e for the electron and j 
= i for the ion), Alfven Mach number {Ma) = 174, magnetosonic Mach number (Mg) = 
159, frequency ratio of the electron plasma oscillation to the electron cyclotron oscillation 



(cjpe/^ce) is 120, and mass ratio M/m = 100. Where fice = eBQ/mc, ujpe = \J ATrriQe^ / m so 
that (jjpe/VLce oc ^Juq/Bq. The quantities wq, e, M, m, and c are respectively, the upstream 
density, the electric charge, the ion and electron mass, and the speed of light. The elec- 
tron plasma oscillation time {2txu^^ = Tpg) is divided into 115 numerical time steps. Each 
simulation cell has 150 particles for each species of the electron and the ion. 

We focus here to follow overall dynamics of shock propagation as well as accurate electron 
dynamics under some realistic plasma parameters, such as high Mach number and the high 
frequency ratio of Wpe/^ce, with statistically enough particles. Under realistic mass ratio, 
we require upstream incoming plasma flow of ~ 6000 km/s for above setup condition of 
ujpe/^ce = 120 and = 174 which is a rather strong shock wave but not unrealistic (for 
example, at RCW 86 shock speed is estimated to 6000 ± 2800 km/s in Ref.|20|). Although, 
our simulation has been carried out under ID and reduced mass ratio condition, results will 
be one of the useful preparation steps toward complete elucidation of the real supernova 
shock dynamics. About some two-dimensional effect for the electron dynamics, although in 
different parameter ranges, there are, for example, Amano and Hoshino |9] for non-relativistic 



case and Kato 



2l|, Spitkovsky [22] for relativistic case. 



III. SIMULATION RESULTS 



A. Overview of shock wave properties 

Figure 1 shows a snapshot of a part of the simulation box in which a shock wave is 
propagating toward the left. From top, the electron momentum phase space density of x 
and y direction {Pf,^ and Pey)-, ion momentum phase space density of x and y direction 
{Pix and Piy)i electric field E^, and magnetic field B^. These quantities are normalized, 
respectively, by muQ for the electron, Muq for the ion, by the upstream motional electric field 
Eq = uqBo/c, and by the upstream magnetic field Bo. The spatial unit for the horizontal axis 
is uq/Qci (= Lstr) which is a characteristic scale of the STR, where VL^ = cBq/Mc is the ion 
cyclotron frequency. Some arrows labeled (a) ~ (e) will be referred to discuss the electron 
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energy spectra in subsection III.D and Figure 8. The profile shows large and intense 
fluctuation generated by microscopic plasma instabilities especially around x = 0.8 ~ 1.2 
(see subsection III.C in detail). 

Since magnetic field profile looks unfamiliar compared with a classical shock profile, 
in Figure 2, smoothing profiles (use right tics) are drawn over raw-data profiles (use left 
tics) with time interval ~ O.lfi"^. The middle panel is a snapshot at the same time of 
Figure 1. The time duration from bottom to a top panel corresponds to about 1.4 Tstr 
(= Lstr/us, characteristic STR passage time), where Us = l.5uo is an average incoming 
flow speed in the shock frame, i.e., shock speed. The normalization scales Lstr and Tstr 
correspond to 115.2c/upi and 8000;"/, respectively, where Upi = ^/^imQe^^jM is the ion 
plasma frequency. Smoothed proflles with low-resolution are obtained by spatial average 
over 5 c/upi with each drawing point shifting at an interval of c/ujpi. Smoothed profiles help 
to identify familiar shock structure with amplified downstream field about three times as 
large as upstream magnetic field as expected from Rankin- Hugoniot (R-H) relation (in our 
ID case, compression ratio is three). As a result of the strong dissipation (thermalization) 
process as shown in the next subsection, in a low-resolution view, the field profiles and 
^article distributions (not shown) show little qualitative variation during their propagation 



23l . |24| |. We can discuss properties appear in Figure 1 (normalized simulation time is 7.1) 



as a typical state of the shock wave. Figures 3-6, 8, and 9 are the same time with Figure 1 
and Figures 7, 10, and 11 include data from t = 7.1 ~ 7.7. 



B. Rapid energy re-distribution 

Rapid and large velocity spread of the electron is seen at just entrance of the STR 
{x = 0.8 ~ 1.0) in Figure 1. In Figure 3 the quantity of the electron thermalization is shown 
as the effective temperature, Te with solid lines compared to the shock flow energy Mu^ 
(top: linear scale) and to the initial temperature T^q = T^q = Tq (bottom: log scale). The 
effective temperature is obtained by calculating f f{p)'ym{V — V)'^dp/ f f{p)dp, where / 
is a distribution function with p = '-fmV, 7 is the Lorentz factor, V is the particle velocity, 
and V is the averaged particle velocity. The ion temperature (Tj) calculated in the same way 
is also plotted with dotted lines in Figure 3. The local enhancement of the Tj at x ~ 0.85 
comes from non-Maxwellian behavior with superposition of the incoming and the reflected 
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ion components. When Tj enhancement appears due to the ion reflection, raise up above 
Tj (> iC'Teo) within a scale of Q.ILstr (leading edge of a classical foot region). In the past, 
a simulation of mildly relativistic low Mach number shocks {Ma ~ 2) under low frequency 
ratio condition (wpe/^ci < 1) shows quick electron energization, which may realize energy 



equipartition between ions and electrons (e.g. Ref.j25|]). To our knowledge, however, such 
strong electron thermalization observed in the present simulation (Tg > Tj) has not attained 
in any other nonrelativistic, lower Mach number shock wave simulations. The remarkable 
abrupt increase of Te indicated by arrow at x ~ 0.92 corresponds to the position of an onset 
of the dynamical change caused by the nonlinear saturation of the plasma instability (see 
next subsection). This means that the nonlinear saturation process plays an important role 
in the rapid energy equipartition between ions and electrons. The local hump of Tg > Tj 
around a; = 1.5 ~ 1.9 is produced intermittently by a temporary variation of the macroscopic 
cross-shock potential structure. In Figure 4, profiles of (middle, normalized by Eq) and 
the cross-shock potential (bottom, normalized by Mus'^/2e) are shown after smoothed over 
ion-scale {25c/upi for E^, 5c/upi for the potential). Overall variation of E^ is brought by 
the polarization between the bulk electron and bulk gyrating ion. The bulk electron is 
overshooting the bulk ion motion. In the top panel of Figure 4, bulk flow variation of the 
ion (solid), the electron (dashed) are shown with the electron excess amplifled 10 times 
(dotted). Several arrows indicate directions overshooting electrons converge to or diverge 
to, which causes Ex variation and resultant cross-shock potential. The cross-shock potential 
is one of the fundamental shock structure 26| and also a critical quantity in terms of the 



electron injection efficiency in DSA process |l4j. l27l|. Time-averaged cross-shock potential is 
about 0.3 ~ 0.4MM,V2e in our case and consistent to the value used in Amano and Hoshino 



Figure 5 illustrates spatial variation of the energy density distribution (% in a logscale) 
among the (a) particles, (b) magnetic field B^, and (c) electric field E^. In the panel (a), 
the ion energy is drawn by dotted line and the electron energy is by solid line. The electron 
energy is calculated by (7 — l)mc^ (for ion energy, m is replaced by M). Each quantity is 
normalized by the total energy of all components. Initially, the energy density ratio is about 
99% for the ion and about 1.0% for the electron. Almost half of the energy initially carried 
by the incoming ion is converted quit rapidly into the electrons within a scale < O.ILstr- 
The maximum energies of and E^ are only of the order ~1% and ~0.1%, respectively. 
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Although, the energy carried by the fields is much less than the particle energy, the fields 
are indispensable especially to the electron energization process as shown below. 



C. Structure of the shock transition region (STR) 

Figure 6 shows an enlargement picture of the entrance of the STR {x = 0.8 ~ 1.15) 
with Pix, E^, the electric current of y-component Jy (in a online color version, the total 
current is black line, the electron current is red line), and B^, from top to bottom. This 
region can be, roughly speaking, divided into two regions by the vertical line at x =0.92. 
Following the change of the ion phase space distribution, other physical quantities also show 
clear qualitative change around this line. The left region (we call here "region I") is where 
we can easily distinguish the incoming ion (positive velocity component) from the refiected 
ion (negative velocity component) in the Pi^ panel. In the region behind the region I (we 
call here "region II"), the mixing between the incoming and the refiected ion population 
progresses rapidly in the phase space. Let us point out some conspicuous structures seen 
in these regions: the leading packet of the magnetic field (region I) and following magnetic 
solitary humps with large amplitude electric field oscillation (region II). In the region I, the 
amplitude of the magnetic field oscillation reaches 10 times as large as the averaged value 
obtained from R-H relation and in the region II, the amplitude of the electric field oscillation 
has over 100 times as large as upstream motional electric field (Eq). In the region I, the 
nonlinear evolution of the two-stream instability between the decelerated hot electron and 



the incoming ion [28| causes some fiuctuation on the ion-acoustic mode. As a result, a series 
of electrostatic potentials are generated which is seen as modulation of the incoming ion as 
well as the electric field oscillation. Some electrons are trapped around separatrix of these 
potentials and, consequently, dragged toward +y direction due to motional electric field. 
Untrapped electrons in non-separatrix region tend to drift toward —y direction under the 
infiuence of local E x B drift. As a result, a series of clear electric current filaments are 
formed by the trapped and untrapped electrons (spiky variation from negative to positive 
in the Jy panel). Since the partial derivatives along y and z vanish in our ID geometry and, 
in the Ampere's circuital law, the Jy term is more dominant here than the displacement 
current term dEy/dt, the relation reduces to —\dB;^/dx\ oc Jy, which facilitates the spatial 
modulation of in the leading packet. A series of resultant negative-positive structures 
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of the magnetic field indicates tlie possibility of the magnetic reconnection when we carry 
out 2D or 3D simulations. The leading packet propagates toward upstream with a velocity 
nearly equal to the shock speed, while each wavy component of a large-amplitude magnetic 
field drifts toward the downstream in the shock frame. 

After the electron is heated rapidly and the ion components begin to be diffusive in the 
region I, another instability between the incoming and reflected ions (ion-ion instability 
hereafter) evolves. In the region II, mixing between the incoming and the reflected ion 
component becomes almost complete. The trapped electrons in the region I are detrapped 
here and clear Jy filaments and negative B-, are disappear in the region II. During this 
process, some components in the leading packet are modulated into the magnetic humps. 
A possible cause of the magnetic hump formation is evolution and saturation of the ion- 
ion instability, not general ion-ion "stream" instability (e.g. Wu et al. 29| in ID, Ohira 



and Takahara 



30| in 2D). This ion-ion instability looks nonlinear direct coupling between 



the incoming and reflected ion modulation which is supported by the hot electron on the 
ion-acoustic mode. The linear analysis shows two-way unstable mode (corresponding to 
the incoming and reflected ion) with the wavelength ~ 2c/ujpi which is a typical scale of 
the magnetic humps with rather large growth rate of 20 ~ 30r2ci (see Appendix). Once the 
magnetic hump structures are formed, they survive over a scale of ~ Lstr in the downstream 
region. Some analytical study about the magnetic humps (magnetosonic solitons), although 
under much smaller amplitude condition, is reported by Pokhotelov et al. 31|. 

The Fix panel in Figure 6 shows a series of ion phase holes are generated by the ion-ion 
instability in the region II. In Figure 7, trajectories of the ion phase holes are shown as the 
propagation of accompanying in a x-time phase space. The two-way complex propagation 
of the ion phase holes affects electron energization process in a stochastic way. These two- 
way streaming of the ion phase holes can survive Landau damping because the phase speed 
of the ion phase holes is much smaller than the thermal speed of the hot electron. 

In spite of such a high Mach number shock, Buneman instability between the incoming 



dose not have strong effect on the 



electron and the reflected ion previously reported (e.g. 
electron energization here in comparison with other instabilities. The Buneman instability, 
however, plays an important role to decelerate electrons and make velocity difference from 
the incoming ion population at the leading edge of the STR. 

When plasma /3 is changed to a larger value, for example, ~ 1.0, we confirm that a 
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structure of the STR is almost identical to ^ ~ 0.1 case here (not shown). The STR structure 
depends strongly on the separation in the velocity phase space between the incoming and 
the reflected ion components and on their thermal width, because that affects the way of 
free energy release via evolution of the plasma kinetic instabilities. At a high-Mach number 
regime treated here, the velocity separation is large enough to trigger strong nonlinear 
evolution of the two-stream instability, whichever j3 may be 0.1 ~ 1. 

D. Electron energization process 

As seen in Figure 5(a), the electron energy density becomes comparable to the ion energy 
density instantaneously. How much energy is attained by the most energetic electron at the 
STR? Figure 8 shows spatial evolution of the electron energy spectra around the STR: at 
just entrance of the STR (a. x = 0.85 ~ 0.95), around the ramp region (b. x = 0.95 ~ 1.1), 
around the magnetic overshoot region (c. a; = 1.1 ~ 1.3), around the magnetic undershoot 
region (d. x = 1.3 ~ 1.6), and its downstream region (e. x = 1.6 ~ 1.9). These regions 
corresponds to the arrows from (a) to (e) in the top of Figure 1. In Figure 8(A) & (C) the 
energy is normalized by the injection ion flow energy in the downstream frame (Mmq/2). 
The panel (B) shows spectra in terms of 7. The spectrum becomes harder as goes deeper 
in the STR. It is surprising that some electrons gain energy more than several times of the 
incoming ion energy just within the entrance of the STR. The most energetic electrons gain 
energy about 1.5 times larger than the kinetic energy gain due to the ion specular reflection 
at the shock front. Each spectrum shows strong thermalization and a little enhancement 
of energetic electron tails. In order to compare thermal population to the Maxwellian, we 
show the energy spectrum at d. with a corresponding Maxwellian of the same effective 
temperature of the region x = 1.3 ~ 1.6 in the panel (C). Although we do not find a 
clear non-thermal tail, some excess is confirmed in higher energy region compared to the 
Maxwellian. 

In Figure 9, the electron number distribution of different energy ranges are shown from 
X — 0.7 to X — 1.9. The top panel shows number of electrons with energy range of 1.5 < 

7 < 1.7 (solid) and of 7 > 2.2 (dashed). The middle panel shows population of 7 < 1.1. For 
reference, the smoothed density profile (normalized by the upstream value) is drawn in the 
bottom panel. Note that y-axis of the top panel is a log-scale. The bulk electron distribution 
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of 7 <1-1 dose not entirely follow the shock-compressed density profile. To keep the shock- 
compressed density, more energetic electrons are needed to appear behind x ~ 0.9. The 
lowest energy electrons (7 < 1.1) shows a little but clear increase with a peak at the leading 
packet. More energetic electrons in the top panel show exponential rise- up and keep on their 
flux in the downstream side. The width of the rise-up scale is comparable to the Larmor 
radius (r^) of the energetic electrons [tl — 0.2 for 7 = 2.2 electrons, tl — 0.11 ~ 0.14 for 
7 — 1.5 ~ 1.7 electrons). The energetic particles are confined tightly in the STR because 
the magnetic field has the perpendicular geometry to the shock normal direction. In our 
simulation, contrary to expectation from DSA mechanism, the energetic electrons are not 
hovering around over the upstream and downstream regions. 

Next, electron trajectories arc investigated. Almost all particles experience both of ac- 
celeration and deceleration in a long time duration. We focus here acceleration pattern and 
pick up two trajectories typical in the leading edge of the STR (Figure 10) and its down- 
stream region (Figure 11). In the right panel of Figure 10, the electron trajectory is drawn 
over time-stacked magnetic field {Bz) profile. The strength of Bz is showed by contour bar. 
A background thermal electron injected into the front of the STR gains energy due to the 
trapping by the leading packet [t = 0.0 ~0.066) as well as due to bouncing motion between 
the leading packet and magnetic hump {t > 0.066) up to 7 ~ 3 (about 200 times as large 
as the initial energy). The kinetic energy gain is brought non-adiabatically. The leading 
packet is propagating toward upstream as a whole structure (an array of stripes from the 
lower-right to the upper- left in Figure 10), but each large-amplitude wavy structure is prop- 
agating toward downstream (each stripe is toward the upper-right). On the other hand, a 
series of magnetic humps are almost standing in the electron bulk flow. As a result, the 
large- amplitude wavy structure is propagating toward the humps. A electron is confined and 
gains energy through multiple nonlinear interaction with these converging large amplitude 
Bz structure similar to the Ist-order Fermi type acceleration mechanism. 

Figure 11 shows another electron orbit among a series of the magnetic humps in the same 
format with Figure 10. The energization process here is rather mild and stochastic. We 
confirm that the kinetic energy gain is also brought in a non-adiabatic manner. A rather 
energetic electron has a large Larmor radius, the effect of the fields can not be symmetric for 
a round-trip of the electron gyro- motion so that stochastic energy change is brought about. 
In Figure 11, the electron shows unmagnetized behavior around time t = 0.05~0.06 and 
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t = 0.1~0.15 (indicated by arrows in the right panel) corresponding to the turning point 
of the gyro-motion. The gyro-orbit in a ~ field has a quite small curvature so that 
the electron is going along y-direction for a while. The magnetic mirror force of the humps 
also keeps the electron in the ^ region in this case. After detrappcd from 5^ ~ 
region and magnetized, the electron gains some energy due to the motional electric field like 
so-called "pick-up ion" process. After that, the electron is bouncing between humps and 
gain energy gradually up to 7 ~4. The nonlinear interaction with these standing magnetic 
humps do not, however, always bring energy gain. It depends not only on the local variation 
of the electric fields (see the middle panel of Figure 4) but also dominantly on the 
accumulation of the local Ey variation caused by the bulk flow and variation. Around 
the region shown in Figure 11, there is still positive and negative bulk flow (= u) variation 
due to the macroscopic ion gyro-motion as shown in the top of Figure 4. The variation of 
Bz (seen as amplitude variation of the magnetic humps) is connected with the bulk flow 
variation. The gyro-motion of the energetic electrons can be synchronous with the motional 
electric field Ey — —e{±u)Bz along y-direction. As a result, there is sometimes net energy 
gain through a round-trip of the gyration. 



IV. SUMMARY & DISCUSSION 



The numerical simulation of a coUisionless shock wave is reported under high Mach num- 
ber and high frequency ratio condition {Ma — 174 and Upe/Qce — 120) similar to shock 
waves generated by the strong blast wave of supernovae. We focus on the electron dynamics 
through the nonlinear particle-field coupling process. We found that the first rapid electron 
acceleration occurs in the thin shock front region. Nonlinear evolution of the plasma in- 
stabilities between the ions and electrons and following the ion-ion instability causes strong 
particle energization. At the moment of the nonlinear saturation of the instability, strong 
mixing occurs between the incoming and refiected ion component in the velocity phase space 
as well as some considerable change in the field structures, for example, from the magnetic 
leading packet to the magnetic solitary hump. The leading packet consists of a series of 
large-amplitude negative-positive magnetic fields. (When a 2D or 3D simulation is carried 
out, wc can expect occurrence of the magnetic reconnection there. It may bring about more 
interesting electron-field dynamics.) Since the each component of the leading packet and the 



12 



magnetic hump are converging, some electrons scattered and bouncing between these two 
structures gain energy rapidly by the Ist-order Fermi type process. Some electrons around 
magnetic humps also gain energy in a rather stochastic manner. In the magnetic hump 
region, there is still non-zero bulk flow variation due to the ion gyro-motion. The electrons 
gain energy not only due to the local E^j. resulted from nonlinear instability evolution, but 
also due to the motional electric field Ey accompanied by the magnetic hump Bz and bulk 
flow. Depending on the gyro-phase, an energetic electron sometimes gains energy when its 
gyro-motion couples with the Ey variation of the ion-scale. 

It is surprising that the electron energy density arises rapidly at just an entrance of the 
STR and becomes comparable to the ion energy density. This energy equipartition is pushed 
by the nonlinear saturation process of the plasma instability. So-called "injection problem" in 
DSA mechanism may be settled if our result remain still true under the realistic dimension 
and mass ratio condition. We showed that in our simulation, the energetic electrons are 
confined within a STR while they are interacting with the fields and gaining energy. So 
that we can estimate the maximum energy jmax by equating the electron Larmor radius and 
characteristic scale of the STR, 

mc^ 1—2 — — r iro% . s 

eB ^^"""^ ^ 2 ^' 

where, Fq is incoming fiow Lorentz factor. Using mq = 2/3VaM4, above equation results. 



When we adopt above equation to the shock wave case in this paper, we have 7max ~5 
which is consistent to our result of 7max ~4. With realistic mass ratio Mjm = 1836, 'jmax 
can reach up to 21 (corresponding to ~ 10 MeV). Under a larger mass ratio condition, the 
STR becomes wider in terms of the scale of electron dynamics and the ion inertia becomes 
relatively larger. Since the ion inertia plays an important role to amplify the fields discussed 
in the current paper, we can expect more electron energization. In fact, some periodic 
simulations in Ref. 28|] shows that when the mass ratio is larger up to the realistic mass 
ratio, the electric field supported by the ion inertia becomes more important in the electron 
energization process. 



^lecent observation shows shocked temperature Tg < Tj for strong shocks (e.g. Ref. 
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36|). Contrary to the observation, our simulation result shows Tg > Tj. One reason 
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for that is a scale difference. Our simulation scale is only about several Lstr much smaller 
than observation can resolve, where Lstr = u-o/Qci ~ 4.6 x 10^ km (for interstellar density 
of 0.1 cm~'^), ~ 1.5 X 10^ km (for 1.0 cm~'^) using Upe/Qci = 120 condition. The temperature 
far outside of the STR is, unfortunately, beyond a scope of the paper. Another reason is a 
shock configuration treated here, namely, the ID and perpendicular shock system which is a 
rather strong restriction. For example, ID simulation tends to overheat electrons compared 
to the two-dimensional simulations |9|, l37|. Although strong thermalization indeed occurs in 
our simulation, the hot electrons shows non-Maxwellian distribution and non-thermal part 
grows as leaving off the STR. The hot electrons can be not only a seed of the future non- 
thermal population but also contribute to further acceleration by keeping the field structures 
within a scale ~ Lstr- 

Unlike the classical diffusive shock acceleration where the electrons are hovering around 
shocked and unshocked region due to scattering by large scale MHD waves, the electrons in 
our simulation experience multiple nonlinear interaction, even within one gyro-motion, with 
the large-amplitude field structures which generated through the nonlinear evolution of the 
plasma instabilities. These nonlinear transportations may lead a breakout of the classical 
diffusion model. Recently, without the classical diffusion model, Malkov and Diamond js^ 
also discuss particle acceleration process under multiple nonlinear interaction among a train 
of magnetic structure "shocklets" of c/upi scale in the confined region of the STR (shock 
precursor). As a generator of the structure "shocklets", nonlinear evolution of the plasma 
instability is also considered between cosmic ray flow and shock incident flow. Not only for 
our rather lower energetic regime, also for higher energetic regime, such microscopic coherent 
structures due to the nonlinear evolution of the plasma instabilities confined in the STR is 
expected to play a key role in the particle acceleration process. 
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Appendix A: Linear analysis of the plasma instability 



In order to find finear growth modes of the time scale larger than the electron gyration 
time, we examine the standard plasma dispersion relation under the magnetized electron and 
unmagnetized ion condition. As seen in the shock transition region, two ion components of 
the incoming ions and reflected ions are adopted. We assume that 45% of the incoming ion is 
reflected at the shock front (the reflection ratio is 0.45). The relative drift velocity between 
the two ion components is assumed as 0.5mo(= ^od)- The net electric current is canceled by 
the total flows of the two ions and electron components. As an initial temperature condition, 
we take IxIO^Tq for the incoming ions and SxIO^Tq for the reflected ions. The magnetic 
held is assumed to be ZBq. 

When considering waves perpendicular to the magnetic fleld, the plasma dielectric tensor 
£(fc,a;) is reduced to have the form, 

up' ^ up up 

Sxx = 1 + /„exp(-?7)— h ^CiO^(Cio) + ^Or^(Cir) 

^ UJ^ T] mlce — OU UJ'' UJ^ 

. , , a; 



^yy = ^ + ^^Inew{-v) 



up nD,ce — ^ 

'^■nc. .X , X UJ 



u;^ m Ice — ^ 

where, /c, uj are the wave number and wave frequency with the Doppler shift included. The 
wave frequency in the rest frame {lo^) is deflned by 0;^ = u; + A;Mod(l — 20;^)- ^pei ^pio, ^pir 
are the plasma frequencies for the electron, incoming ion, and reflected ion, respectively. /„ 
and /„' are the modifled Bessel function of n-th order with the argument rj — v'^k'^/2Q'^ 
and its derivative. We take the summation for n from —10 to 10 here. Z is the plasma 
dispersion function with the argument Qq = {u — kuQ^) / kvio for the incoming ion component 
and Qr = (cu + kuQd)/kVir for the reflected ion component. Where ViQ, fj,., and are, 
respectively, the incoming, reflected, and electron thermal velocity. Using above dielectric 
tensor, the dispersion relation is 

2" 



^yy 1 



rck 



+ < = (Al) 



Figure 12 shows some of the solutions of eq.(Al). Top panels include real part of u 
and bottom panels include imaginary part of uu (= 7, growth rate) versus wavelength A. 
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The panels (a) and (b) belongs to the incoming ion because their phase velocity is the 
same direction with the incoming ion. The panels (c) and (d) belongs to the reflected ion 
with the same reason. Their harmonics in other wave number range is not shown, uj and 
A are normalized by the plasma frequency of the upstream ion cjpj and its inertia c/wpj, 
respectively. The electron temperature Te is set up to 2xlO'^To (dotted), 4xlO^To (long- 
dotted), SxlO^To (dashed), and IxIO^Tq (solid). As the electron temperature rises up, the 
growth rate increases with its peak moves a little toward smaller A. The solution does not 
have a strong dependence on the ion temperature but smaller ion temperature results larger 
growth rate. The wavelength ~ 2c/upi (comparable to the normalized width Ax ~ 0.01) 
around the peak of the growth rate corresponds to the width of the magnetic hump stated in 
subsection III.C. Coupling between such linear growth modes in the incoming and reflected 
ion component can be some cause of the magnetic hump generation. However, since strong 
nonlinearity is observed in the simulation, especially in region II, we can not separate clearly 
the linear growth from other nonlinear plasma actions. 
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0.5 1 1.5 2 2.5 



FIG. 1: (Color online) The structure of a high Mach number shock wave. From top, the electron 
momentum of x and y components {Pex and Pey), the ion momentum of x and y components {Pix 
and Piy), the electric field (Ex), and the magnetic field (Bz). The contour shows 10 logarithmic 
distribution of the particle number. The electron and ion momenta are normalized by muo and 
Muq, respectively, the electric and magnetic fields are normalized by the upstream motional electric 
field Eq = uqBq/c, and upstream magnetic field Bq, respectively. The horizontal axis is normalized 
by uq/^Ici- Several arrows (a)-(e) are referred in Figure 8. The normalized simulation time for the 
figure is 7.1. 
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FIG. 2: Bz profiles with raw value (left tics) and with smoothing value (right tics) are drawn at 
an interval of ~ 0.10"^ from bottom to a top panel. Smoothed profiles are obtained by spatial 
average over bc/upi with each drawing point shifted at an interval of c/upi. The time of the middle 
panel corresponds to the time of Figure 1. 
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FIG. 3: The effective temperature variation normalized by the shocked flow energy Muj (top) and 
normalized by the initial temperature Tq (bottom), where Ug (= 1.5uo) is an averaged injection 
flow speed in the shock frame. The dotted and solid lines correspond to the ion and electron 
temperature, respectively. The arrow in the top panel at a; ~ 0.92 shows a position of abrupt rise 
of the Tg. This position corresponds to the boundary between region I and II discussed in the 
subsection III.C. The time of the figure is same with Figure 1. 
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FIG. 4: Top: The bulk flow speed of the ion (solid), electron (dashed), and the flow speed difl'erence 
of these two components (dotted, amplified 10 times) normalized by uq. Arrows in the panel shows 
the direction of the electron component shooting out of the ion component, middle: Ex profile 
normalized by Eg. Top and middle profiles are smoothed over 2bc/upi with each smoothing point 
at an interval of 2c/u)pi. bottom: The electrostatic potential corresponding to the above E^ profile. 
The data is normalized by Mu1/2e and smoothed over bc/ujpi with each smoothing point at an 
interval of c/wpj. The time of the figure is same with Figure 1. 
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FIG. 5: The variation of the energy ratio (%) among (a) the ions and electrons (dotted and sohd 
hne), (b) the magnetic field , and (c) the electric field from the upstream to the downstream side 
of the shock transition region (x = 0.5 ~ 2.5). Note the vertical axis is a log-scale. The time of 
the figure is same with Figure 1. <pq 
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FIG. 6: (Color online) The enlargement picture of Figure 1, an entrance of the shock transition 
region. From top, ion phase space density (Pix), Ex, electric current Jy of the electron (the lower 
line) and total (the upper line), and Bz- Smoothed profiles {Jy and Ex) are obtained by averaged 
over c/ujpe with each drawing point shifted at an interval of c/Qupe- The vertical line indicates 
boundary between the region I and II. The time of the figure is same with Figure 1. 
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FIG. 7: (Color online) Trajectories of large-amplitude electric field {E^) in the x-time phase space. 
The relative amplitude of is shown by the contour bar. Time is normalized by Tstr- The origin 
of the time in the figure is same with Figure 1. 
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FIG. 8: (Color online) The electron energy spectra. The energy is normalized by the ion injection 
flow energy (Mmq/2) in the panels (A) & (C). Each spectrum in the panels (A) & (B) is sampled, 
from left, at just entrance of the STR (a. x = 0.85 ~ 0.95), around the ramp region (b. x = 0.95 ~ 
1.1), around the magnetic overshoot region (c. x = 1.1 ~ 1.3), around the magnetic undershoot 
region (d. x = 1.3 ~ 1-6), and its downstream region (e. x = 1.6 ~ 1.9). These regions corresponds 
to the regions indicated by arrows (a) ~ (e) in Figure 1. In the panel (C) the energy spectrum c. 
is drawn with a corresponding Maxwellian with the same effective temperature in this region. The 
time of the figure is same with Figure 1. 




FIG. 9: The electron number distribution with energy 7 > 2.2 (top dashed Une), 1.5 < 7 < 1.7 
(top sohd Une), and 7 < 1.1 (middle solid line). The electron density profile normalized by the 

upstream density is shown for reference (bottom). The density profile is smoothed by averaged 
over c/ujpi. Note the top vertical axis is a log-scale. The time of the figure is same with Figure 1. 
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FIG. 10: (Color online) An example of the electron orbit around the entrance of the STR on the 
time-stacked profiles of the magnetic field Bz (right). The strength of Bz is shown by the contour 
The left panel shows the electron energy variation. The origin of the time in the figure is same 
with Figure 1. 
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FIG. 11: (Color online) An example of the electron orbit and energy variation in the downstream 
side. The format is the same with Figure 10. Two arrows in the right panel indicate the regions 
where the electron shows unmagnetized behavior. The origin of the time in the figure is same with 
Figure 1. 
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FIG. 12: The dispersion relation oj versus A (wavenumber) ((a) &; (c)) and growth rate 7 versus A 
((b) & (d)) of the peculiar mode in the regime of the ion- ion interaction on the incoming ion ((a) 
& (b)), and on the reflected ion ((c) k. (d)). uj and A are normalized by the plasma frequency of 
the upstream ion ujpi and its inertia c/upi, respectively. 
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